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I. Introduction 


This is the final contract and scientific mission analysis report 
required under contract number NAS 9-10579 between the NASA and Leland 
Stanford Junior University. This work was suggested in a proposal that 
S-band (13 cm) radio frequency transmissions from the lunar orbiting 
Command Service Module be used to obtain bistatic radar echoes from the 
lunar surface* The reception of these echoes was to be carried out on 
the earth. The proposal was accepted and a contract initiated with an 
authorization to incur cost dated January 29 , 1970- Tiie P er i°d of per- 
formance was from that date through June 30 , 1973* 

Subsequently 9 it was realized that VHF (ll 6 cm) radio frequency 
transmissions were also available and could be used for bistatic radar 
purposes* Due 10 the straight forward nature of the experiment and the 
fact that no special on-board equipment was required, it was agreed that 
the bistatic experiments would be carried out as early as possible in 
the mission set. The initial experiment schedule included VHF observa- 
tions on Apollo 13 , and both S-band and VHF observations simultaneously 
on Apollos lk and 15 . These observations were later expanded to include 
the simultaneous dual wavelength experiment on Apollo 1 6 as well* No 
observations were obtained from Apollo 13 due to the unfortunate circum- 
stances associated with that mission. Simultaneous S-band and VHF 
observations were obtained on each of the three subsequent missions of 
Apollos 14, 15 and 1 6 . 
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I X . Experiment Purpose 


The purpose of S-I 7 O was to: 

1 # Determine the principal electromagnetic parameters of the lunar 
crust ; 

2 Determine the electromagnetic properties of lunar regolith and 
subsurface structure to depths of at least several tens of 
meters ; 

3 . Determine the frequency dependence, if any, of the principal 
crustal electromagnetic parameters ; 

4. Determine the presence of vertical surface structure and spatial 
distribution of horizontal surface structure within the range 
given above with emphasis placed on the depth of the regolith. 

5 . Map the statistical properties of surface roughness on the 
scales of the wavelength employed. 

In general, we have attempted to develop quantitative measures of lunar 
surface properties that can be used to characterize and define lunar 
surface structure. These results have then been applied to the larger 
problem of lunar evolution and state through collaboration with colleagues 
in the geological and geophysical sciences, 

III. Publications 

To date this experiment has resulted in the following publications: 

Tyler, G. L. with Thompson, T. W. et al., "Remote Sensing of Mare 

Serenitatis, " to appear in Apollo 17 Preliminary Science Report , NASA. 

Tyler, G. L. with Moore, H. J. "Comparison Between Photogrammetric and 
BistatiC“Radar Slope-Frequency Distributions," to appear in Apollo 17 
Preliminary Science Report, NASA. 

Tyler, G. L. and Howard, H. T., "Dual Frequency Bistatic Radar Investiga- 
tions of the Moon with Apollos 14 and 15 , " J. Geophys. Res ., 
to appear August, 1973* 

Parker, M. N. and Tyler, G. L., "Bistatic Radar Estimation of Surface 
Slope Probability Distributions with Applications to the Moon," 

Radio Science, 8 , PP . 177-184, March, I 973 . 


'1 


2 



Howard, H. T. and Tyler, G. L., "Bistatic Radar Investigation," Apollo 16 
Preliminary Science Report , NASA Publication SP- 315 , pp. 25 - 1 _ - 25-13,* 

1973 . 

Howard, H. T. and Tyler, G. L., "Bistatic Radar Investigation," Apollo 15 
Preliminary Science Report, NASA Publication SP- 289 , pp. 23 - 1 23-10, 
1972. 

Howard, H. T. and Tyler, G. L., "Bistatic Radar Studies of the Lunar 
Surface," Apollo lU Preliminary Science Report. NASA Publication 
SP-272, pp. 247-266, 1971. 

IV. Experiment Discussion and Results 

Observational data records, reduced data, and a detailed report of 
the data reduction process from this experiment have been submitted to 
the National Space Science Data Center. A copy of the NSSDC report 
(Stanford-Apollo Bistatic-Radar Experiment (S- 170 ) : National Space 

Science Data Center Data Description) is bound with this report. It 
contains the principal technical details of the experiment. 

Additional observations were carried out using the ll 6 cm wave- 
length alone. Approximately three times more data were obtained in this 
mode than in the dual wavelength case. The single wavelength observa- 
tions represent an experiment of opportunity in that by agreement with 
JSC the experimenters exercised no control over the lunar ground track, 
the time of observation, or the spacecraft antenna directivity with 
respect to the experimental geometry. Consequently, these observations 
are somewhat less useful than those obtained at the dual wavelength in 
terms of direct in inference of surface properties. However, since data 
were obtained on sequential orbits, a two dimensional scanning action 
across the lunar surface was obtained. Thus, these data are quite 
valuable for mapping two dimensional variations in .lunar scattering 
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propel ties . Due to the heavy emphasis on the dual wavelength observa— 
tionSj these data have received little attention thus far. 

This experiment grew from earlier work with Lunar Orbiters I, III, 
(Tyler et al. 1967) and Lunar Explorer 35 (Tyler 1968) . The Lunar 
Orb iter observations provided the first bistatic-radar echoes from the 
m °on> demonstrated the feasibility of such experiments, and indicated 
that oblique incidence forward scattering from the lunar surface was 
dominated by quasi-specular scattering mechanisms. The first systematic, 
scientific investigation of an extraterrestrial body using the bistatic 
radar technique was carried out with Explorer 35 spacecraft. Only a 
single wavelength (2.2 meters) was available for these observations. 

The results of these experiments have been reported elsewhere (e.g., 
v. Tyler 1968: Tyler and Simpson 1070 ; Tyler, Simpson and Moore 1971) 

The Apollo observations continued this work. Although the data set 
obtained is more limited with respect to coverage than that from Explorer 
35, the experimental controls were much more stringent and the experi- 
mental signal -to-noise ratio was approximately an order of magnitude 
higher. The simultaneous, dual-wavelength nature of this experiment 
is unique. First order analysis and interpretation of the dual-wavelength 
data have been completed. From the conclusions of Tyler and Howard , 

( 1973 ): 

"The oblique geometry scattering properties of the lunar surface 
are wavelength dependent in the decimeter to meter wavelength range. 
At a fixed wavelength, the scattering law is highly dependent on 
local topography. There are systematic differences in the average 
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scattering properties of mare and highland units . At 13 cm wave- 
length mare surfaces are well modeled by a uniform dielectric with 
relative dielectric constant e = 3-1- Deviations from the value 
at this wavelength with the units studied thus far are small, on 
the order of 10$, At 116 cm the mare have considerable variation 
in inferred reflectivity. While much of the surface can be well 
matched by uniform dielectric models similar to that at 13 cm 
wavelength, there are also considerable deviations from that model. 
No simple model, consisting of a uniform dielectric interface or a 
contiguous set of such interfaces will match the' data. Some type 
of frequency dispersive mechanism, such as layering and lossy 
dielectric is required. Highlands units in the Palus Sornni, and 
Mare Crisium rim area possess lower effective dielectric constant 
than the mare. At 116 cm wavelength the data in this area are well 
modeled by uniform dielectric interface with relative dielectric 
constant e ~ 2.8. At 13 cm wavelength it has not been possible 
to find the simple model for the interface in this area. In the 
Apennines, at approximately 25° incidence, both 13 and ll6 cm data 
show a reduction in reflectivity consistant with a change in rela- 
tive dielectric constant from e = 3,1 to e = 2.8. Similar 
results are obtained in the central highlands. 

Lunar rms slopes and slope distributions inferred from the radar 
data also show considerable variation with probing wavelength and topo- 
graphy. The systematic differences between mare and highlands are 
even stronger. At 13 cm, highland rms slopes vary approximately between 
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tan 6° and tan 8°, while in the ir.are this range of variation is approxi- 
mately tan 2° to tan 4°. At ll6 cm^ highland rms slopes nearly equal 

those of the shorter wavelength, but within the mare the ll6 cm slopes 

tend to be very nearly one-half as large as the 13 cm values at the same 
location. Evidently the small scale roughness is greater in the mare 
than that in the highlands * Changes in the rms slopes within the mare 

are a sensitive indicator of changes in terrain types. In some cases, 

such changes may be interpreted as occurring at flow contacts, or 
between geological units of different age or history. 

In forward scatter, the unpolarized part of the 13 cm radar echo 
is generally insensitive to lunar topography. Variations in the unpola- 
rized part are small and apparently without the systematic behavior among 
units evidenced by the polarized part. At ll6 cm, there are significant 
variations in the unpolarized power and scattering law, but these varia- 
tions are not associated with large scale lunar topography in an obvious 
manner. Chance observations of several craters that are prominent back- 
scatterers in this wavelength range indicate that conversion of polariza- 
tion state is primarily responsible for enhancement of the depolarized 
backscatter radar echoes...." 

The data have demonstrated: 

1. A strong wavelength dependence in the oblique scattering proper- 
ties of the lunar surface in the decimeter to meter wavelength 
range ; 

2. Systematic differences in average scattering properties at a 
single wavelength of mare and highland units. 

3. Variations in lunar surface structure with depth within mare 
units; 
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U. Systematic differences in lunar surface shape, as expressed by 
variation in lunar rms slopes with wavelength and between mare 
and highland units; 

5 . The presence at 13 cm wavelength of a diffuse scattering 
component that is essentially constant over the lunar surface; 

6. The presence at 116 cm wavelength of a highly variable diffuse 
scattering component that is not correlated with obvious sur- 
face units. 

These results have been given quantitative form and are described in the 
Preliminary Science Reports and in Tyler and Howard (1973) • believe 

that while significant results have been obtained, considerable work 
still remains. 

In addition to the first order analysis techniques on which the 
foregoing conclusions were based, we have also developed techniques for 
extracting detailed slope frequency distributions from bi static radar 
data (Parker and Tyler 1973). Preliminary results from this type of 
analysis are given in the paper just referenced, and in the Apollo 17 
Preliminary Science Report (Moore and Tyler 1973)* This technique 
provides a powerful method for examining subtle variations in the geo- 
logic texture of the lunar surface. A further outgrowth of this work 
is the ability to derive quantitative measures of the three dimensional 
distribution of scattered power in a particular scattering area. 

We have also undertaken joint efforts with other investigators, 
principally Henry Moore of the USGS, Menlo Park, and T. W. Thompson of 
the Jet Propulsion Laboratory, Pasadena. An object of this collabora- 
tion was to carry out detailed studies of certain lunar surface features, 
for example Mare Serenitatus, using available data from a broad range 
of remote sensors and Apollo observations. Toward this end, a common 
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data file has been created at Stanford and the USGS in Menlo Park. At 
present, this file contains the dual frequency bistatic radar observa- 
tions and an independently derived geologic interpretation of the lurrain 
for each bi static radar data point. We intend to incorporate other 
observations, such as infrared eclipse temperatures, terrestrial back- 
scatter radar maps, and lunar color differences, into this data file. 

Beyond these preliminary results the Apollo data have provided the 
basis for a considerable advance in radar studies of the moon. The 
lunar surface is well resolved, the data are of high signal-to-noise 
ratio, and the theory of forward scatter from rough surfaces is 
sufficiently well understood to permit detailed inferences regarding 
the properties of the lunar surface in the scattering region. Except 
near the sub-earth point such inferences cannot be drawn from earth 
based radar observations because of the much more complex nature of 
the radar backscattering process at oblique incidence. Lunar rms slopes, 
slope frequency distributions, effective surface reflectivity, and the 
wavelength dependence of these quantities have been obtained from the 
Apollo data. But further work is required to bring these results; into 
focus with the lunar problem as a whole. 

As mentioned previously, the analyses carried out thus far have 
dealt with the global nature of the bistatic radar data. That is, the 
data have been analyzed in terms of their average variations with earth, 
moon, and spacecraft geometry. Large scale lunar features such as mare 
and terra units have been taken into account in the analysis and inter- 
pretation, but the primary aim of the work has been to describe the 
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properties of the moon as a whole. The techniques used to date are also 
suitable for studies of more detailed lunar surface features. In addition, 
there are certain aspects of the data that are not completely understood 
but where we can expect progress. We have proposed that this work be 
continued under the Lunar Analysis and Synthesis Program. 

V. Conclusion 

While we have concentrated here on the scientific results of S-l^O^ 
we wish to emphasize that most of our efforts over the period of the 
contract have been involved in the technical process of collecting and 
reducing the data. These efforts are described in considerable detail 
in the NSSDC Report submitted with this report. The interested reader 
is referred to it for information regarding the nature of these efforts. 

Over Lhe period of this contract we were privileged to work with a 
large number of individuals at JSC and within other NASA organizations. 

We wish to comment positively here on the generally high caliber of the 
personnel with whom we worked. We wish to make special note of the 
helpfulness and cooperation of three individuals, Louis Leopold of the 
Electromagnetic Systems Branch, Patrick Lafferty of the Science and 
Applications Directorate and Harley Weyer of Flight Operations. Their 
roles were vital to the success of this experiment. 

In general, we feel that this experiment has been highly successful* 

It has resulted in a distinct contribution to the radar astronomy art and 
a unique lunar data set. We have established certain results regarding 
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the nature of the lunar surface as a whole . Our remaining tasks are to 
extend this work to the study of specific lunar features, in collabora- 
tion with our geological colleagues, and to apply the technique to other 
celestial bodies. 
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Introduction 


The Stanford Apollo Bistatic-Radar Experiment (S- 170 ) was carried 
out during the lunar orbit phase of the Apollo lh, 15, and 1 6 flights* 

The experiment was unique in that no special equipment was placed 
onboard the Apollo vehicles and carried to the moon, and that the 
principal observations were carried out on the ground. Radio-frequency 
transmissions from the orbiting command-service-module were directed 
toward the moon and received on the earth after reflection from the 
lunar surface. Two wavelengths, 13 cm (2287.5 MHz, S-band) and ll6 cm 
( 259.7 MHz, VHF), were used. The best data were obtained during periods 
when the spacecraft was maneuvered to maintain a predetermined, although 
changing, attitude with respect to the earth and moon. During these 
periods, data were obtained at the two wavelengths simultaneously. 

Data were also obtained at the ll 6 cm wavelength during periods of 
inertial hold and SIM bay attitude maneuvers. Data reception took 
place at two sites on the earth, the NASA-DSN 6 b m antenna facility 
located near Barstow, California and the Stanford Research Institute/ 
Stanford University U 6 m antenna facility located on the campus of 
Stanford University, Stanford, California. The NASA facility was used 
for reception of the 13 cm data, the Stanford facility for the ll 6 cm 
data. In both cases the elements of the receiving system critical to 
the experimental objectives were under direct control of the experimenters . 
At the DSN, a Signal Conditioning Unit designed and constructed at Stan- 
ford was inserted as a critical series element in the data receiving 
system and used to set system bandwidths, levels and timing information 
for the analog recording system. Input signal levels, bandwidths, and 
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frequencies to the Signal Conditioning Unit from the DSN were specified 
and monitored in real-time during data collection by the experimenters. 
Data reduction was carried out at the Stanford Sigma 5 Real-Time Computa- 
tion Facility. With two exceptions (cf. Notes on Block Diagram I, 4 d, 
and Appendix V) identical programs and procedures were used for 13 cm 
and llo cm data. Also, with one exception, the same programs were used 
for all three flights (cf. Subtask 2). However, certain parameters, de- 
fined elsewhere in the report, were varied for data at the two wavelengths. 

Good data were obtained from the three flights. In the data sets 
accompanying this report we include the simultaneous 13 cm and ll6 cm 
observations, but not the ll6 cm data obtained in the inertial hold and 
SIM bay attitudes. Data are given in two forms: 

a) A complete set of observations reduced to short time averages 
of the electromagnetic wave spectra for the 13 cm and ll6 cm 
observations. These observations have been corrected for 
instrumental effects, and are unedited. Tapes containing 
these data are referred to as JM Doptrack tapes. The 
observations have been merged with trajectory data obtained 
from the Manned Space Flight Center in Houston, and certain 
ancillary data computed from the trajectory. 

b) A complete set of reduced data records, called Integral tapes, 
obtained from the JM Doptrack tapes, which describe certain 
properties of the JM Doptrack data, such as moments of the 
echo spectra, and inferred properties of the lunar surface, 
such as rms slopes. 

There is one-to-one correspondence between the Integral data and the JM 
Doptrack data; an interested investigator may use either the reduced 
data records given on the Integral tapes, or use the JM Doptrack tapes 
to verify the reduction procedure and perform additional reductions. 

The remainder of this report describes the data collection and 
reduction system in some detail, attempts to give cautioning notes to 

X5 
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the data user, and describes the tape contents and formats for JM Dop- 
track and Integral tapes. However, no particular attempt to explain the 
motivation for the overall data reduction procedures will be made, except 
as it affects individual steps not described elsewhere. A more general 
understanding of the experiment, its goals, limitations, and preliminary 
results may be obtained from the references listed below: 

Howard, H. T. and G, L. Tyler, ,J Bistatic-Radar Studies of the 
Lunar Surface,” Apollo l4 Preliminary Science Report , NASA 
publication SP-272, 1971 . 

Howard, H. T, and G. L. Tyler, 1r Bistat ic-Radar Investigation,” 
Apollo 15 Preliminary Science Report, NASA publication SP- 289 , 

P. 23-1, 1972. 

Howard, H. T. and G. L. Tyler, Apollo 1 6 Preliminary Science 
Report, NASA publication SP- 315 , p. 25~1, November, 1972, 

Tyler, G. L. and H. T. Howard, tr Bistatic Radar Observations of 
the Lunar Surface with Apollos 14 and IS,” paper presented at 
Third Lunar Science Conference, Houston, Texas, January, 1972. 

Tyler, G. L. and D. H. H. Ingalls, functional Dependence of 
Bistatic Radar Frequency Spectra on Lunar Scattering Laws,” 

J. Geophys. Res ,, Vol. 76 , No. 20, pp. 4775 -U 785 , July, 1971 . 



Synopsis of Observations 

A synopsis of the observations included with this report is given in 
Table I, "Stanford Apollo Bistatic-Radar Experiment Parameters # " The 
definitions of the columns in Table I are given below: 

SPACECRAFT - Designation of flight number 

ORBIT - NASA-MSC orbit numbers , counted from lunar orbit 
insertion 

WAVELENGTH - Either 1 3 cm or ll6 cm, wavelength of electromag- 
netic radiation from command-servic e-module, 
corresponds to frequencies 2287.5 MHz and 259.7 MHz 
respectively , 

ANTENNA - Command-servic e-module antenna in use. Quantities 
in ( ) are NASA-MSC designation. Antenna patterns 
used in data reduction obtained from NAA (1966a, 

1966b, 1969), and MSC (1967). Apollo 14 experiment 
required compromise in spacecraft attitude to 
accomodate patterns from two fixed antennas, 

3 db BEAMWIDTH - One-half power antenna beamwidths 
POLARIZATION - Of comm and- service-module antenna 


GAIN - Command-service-module antenna gain 
POWER - Transmitted power, from command-service-module 


SYSTEM TEMPERATURE - 


Equivalent temperature of receiving system on the 
earth, looking at cosmic background radiation 


perpendicular to plane of galaxy. 
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TABLE I 





Stanford Apollo Bistatic- 

-Radar Experiment Parameters 



SYSTEM 

spacecraft 

ORBIT 

WAVELENGTH 

ANTENNA 3 db 

BEAMWIDTH 

POLARIZATION 

GAIN^’ 

power y 

TEMPERATURE 

Apollo 

14 

25 

13 cm 

Cavity 
backed 
helix 
(OMNI C) 

'6o° 

Right elliptical 
axial ratio ~7 db 

-1,5 db 

~4 W 

2 T°K ± 3 ° 

Apollo 

i4 

25 

ll6 cm 

Scimitar 
(VHF LEFT) 

a 

Linear, 
maintained 
in plane 
of incidence 

~0 db 

~2.5 w 

~1000°K ± 50° 

Apollo 

15 

28 

13 cm 

Steerable 

crossed 

dipoles 

(high gain, 
wide) 

'4o° 

Right circular 
axial ratio <1,0 db 

-1.5 db 

~4 w 

27°K ± 3° 

Apollo 

15 

28 

ll6 cm 

Scimitar 
(VHF RIGHT) 


Linear, 
varies with 
respect to 
plane of 
incidence 

~0 db 

~2-5 w 

~1000°K ± 50° 

Apol lo 

1 6 

4o 

13 cm 

Steerable 

crossed 

dipoles 

(high gain, 

WIDE) 

'4o° 

Right circular 
axial ratio <1,0 db 

-1,5 db 

~4 w 

27°K ± 3° 

Apollo 

16 

4o 

ll6 cm 

Scimitar 
(VHF LEFT) 

Of 

Linear, 
maint ained 
in plane of 

~0 db 

~2.5 w 

~1000°K ± 50° 


incidence 


a beamwidth of scimitar not defined 
g includes circuit losses, nominal values ± 3 db 

y in carrier signal, nominal values, actual values not measured in flight 
6 varies with orbital position of spacecraft, cold sky values given 
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Data Collection, Processing, Reduction 


Block Diagram I, located at the rear of the report, depicts the flow 
of data through collection, processing, and reduction. Annotations and 
notes (circled numbers) give brief comments or descriptive titles as 
guides for reference. The overall data flow is divided into 11 subtasks, 
indicated by horizontal brackets. The remainder of the report is sub- 
divided according to the subtasks. It is assumed that the reader has 
general familiarity with analog-to-digital and digital techniques for 
data reduction and analysis. Critical data reduction parameters are 
given in Table II and antenna parameters are given in Tables III and IV. 

Each subtask description provides a general explanation of that 
subtask function. An attempt has been made to maintain independent 
description, with minimum reference to other functions. Highly techni- 
cal details, such as magnetic tape formats and coordinate descriptions, 
are reserved for the appendixes. In this way a careful reading of the 
subtask descriptions, in connection with Block Diagram I, should provide 
a good overview of the data reduction process. The appendixes may be 
read at a later time, or for details. 

A THOROUGH STUDY OF BLOCK DIAGRAM I PRIOR TO THE USE OF REPOSITORY 
DATA IS VERY HIGHLY RECOMMENDED. 
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TABLE 11 







STANFORD APOLLO DUAL- FREQUENCY BtSTATIC-RADAR DATA SUMMARY 





APOLLO 

14 

APOLLO 

15 

APOLLO 

16 


116 cm (VHF) 

13 cm ( 3 -BAND) 

116 cm (VHF) 

13 cm (S-BAND) 

1I6 cm (VHF? 

13 cm ( 8 -BAND) 

Julian Ephemeris Day 








(00:00 UT2 Preceding 
Data) 

24409&8.5 

244U64.5 


2441430.5 

Refer an ce Epoch for 







. . . 

Coordinate Systems 

21*1+0952.509 

2441317-752- 


2441317-752 

CALENDAR DATE 

February 

6 , 1971 

August 1, 1971 


April 23, 1972 

Orbit Number (From Lunar 




28 




Orbit Insertion) 


25 



40 | 

Transmitter Frequency 

259.7 MHz 

2267.5 MHz 

259-7 “Hz 

2237-5 MHz 

259-7 MHz 


2267-5 MH* 

VT2 START/STOP TIME 

q6*3T : 17/07 : 30:31+ 

06 : 3 S; 30 / 07 : 32: 28 

01 :23 :00/02:23:58 

01 e16»30/02j 29:00 

01:16*30/02:16:26 

01:17:00/02:27:05 

Spacecraft Antenna 
(MSC Designation) 

VHF LEFT 

OMNI C 

VHF RIGHT 

HIGH GAIN (HIDE) 

VHP LEFT 


HIGH GAIN (WIDE) 

High Gain Antenna 
Pointing Angles 




144.° 




or 




122.° 



302.° 

P 

“ 








Stanford 

DSN/STANFORD 

Stanford 

dsn/stanford 

Stanford 


DSN/STANFORD 

Data Source 

10.003 MHz i.f.rcvr 

Closed Loop revr 

10.003 MBe i.f.rcvr 

Closed Loop revr 

10.003 MHz i 

f . revr 

Closed Loop revr 

Receiver Bandwidth 

-30 kHz 

~20.0 kHz 

■*3-5 kHz 

"20.0 kHz 

~3.p kHz 


"20.0 kHz 

Data Sampling Frequency 

10.0 kHz 

43.0 kHz 

10.0 kHz 

43.0 kHz 

10.0 kHz 


43.0 kHz 

Total Analysis Band- 1 - 
width 

5.0 kHz 

21.5 kHz 

5.0 kHz 

21.5 kHz 

5.0 kHz 


21.5 kHz 

N, Number of Analysis 
Bins 

1024 

1024 

2048 

1024 

2048 


1024 

Analysis Resolution 

9.6 Hz 

42.0 Hz 

4.9 Hz 

42.0 HZ 

4.9 Hz 


42.0 Hz 

Length of Data 
Wi nd ow/I rans f om 

0.1024 sec. 

0.02381395 sec. 

0.2048 £ ec. 

0.62381395 sec. 

0.2048 sec. 


0.02381395 »<*. 

L, Number of Transforms 
Averaged per JM Record 
(see Subtask 4) 

26 

100 

23 

99 

23 


99 

Frame Length 

2.6624 BOO. 

2.351395 

4.7104 ^:ec. 

2.357581 »ec. 

4.7104 sec. 


2.357581 sec- 
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TABLE III 


Stanford Research Institute 46 m Antenna Performance (ll 6 cm) 


Efficiency Feed System Mount 

~3 5$ ±5$ Crossed dipole array el/ az 

axial ratio ~ 1 , db ± 0 9 ^ a 
isolation ~l 6 db i 1.0 

TABLE IV 

NASA DSN 64 m Antenna Performance (13 cm)* 

Aperture Efficiency Feed System 

64 m dia. 58 1° ±4^ Waveguide horn 

axial ratios < 0.8 db ^ 
isolation > 26 dbY 

a estimated 

Private communication, D. Bathker, JPL, 1973. 


Mount 

el/az 


Aperture 
46 m dia. 


Pointing Loss 
<0.1 db 


Pointing Loss 
< 0.03 db 


* DSN, 1972 
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Subtask 1 


Data Collection 

a) ll6 cm* A detailed block diagram of the ll6 cm 

receiver is given in Appendix I* The receiver is a super- 
heterodyne of standard design* An unusual feature of this 
system is the summed second local oscillator signal used 
to produce offset 9.O and 10.003 MHz intermediate frequency 
signals. In operation, the 10.003 MHz channel was tuned 
to the downlink signal carrier and its accompanying echo, 
while the 9*0 MHz channel was tuned to the subcarrier signal 
displaced + 31*6 KHz from the 259*7 MHz carrier. This 
procedure was adopted to provide frequency dispersive 
redundancy against locally generated interference at the 
main carrier frequency. The receiver passband characteris- 
tics are given in Appendix III. However, all data given 
here were obtained through the main 10 MHz channel. 

Receiver outputs were multiplexed with standard frequency 
references and clock signals and recorded on analog tape. 
Parameters of the Stanford Research Institute 46 m dish 
are given in Table III. 

(b) 13 cm. A detailed block diagram of the 64 m 13 cm 

receiving system is given in Appendix II. Again, standard 
superheterodyne techniques are employed. All 13 cm data 
were obtained from a phase-locked loop signal tracking 
system using manual tracking when the direct signal dropped 
below threshold of the phase-locked loop. The receiver 
passband characteristics are given in Appendix III. The 
NASA-DSN station configuration for this experiment is 
given elsewhere ( DSN , 1970, 1971 ) . Again, analog signals 
from the receiver output multiplexed with clock and 
reference frequencies were recorded. Parameters of the 
NASA DSN 64 m dish are given in Table IV, 
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Subtask 1 (eont.) 


(c) Critical Bandwidths. Appendix III gives the power spectral 
densities obtained at the receiver outputs for the three 
Apollo experiments with uniform power spectrum (white) 
noise input. These curves were obtained in Subtask 4, 
and were used in the data normalization „ For the ll6 cm 
system additional checks made with coherent signals within 
+ 5 MHz either side of the 10 MHz intermediate frequency 
and within + 20 MHz of the first intermediate frequency 
verified the absence of spurious responses that would not 
be detected by the noise calibration technique. Similar 
tests have been carried out by the DSN and the experimenters 
at the 64 m facility. Thus, the curves in Appendix III 
accurately represent the receiver response to signals. near 
the frequency to which the receivers were tuned. Image 
rejection in the ll6 cm system was greater than 100 db. 
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Subtask 2 


Sampling 

Data sampling was carried out using standard techniques 
on the Stanford Sigma 5 Real-Time Computation Facility. 

Right circular and left circular polarization signal 
channels were sampled simultaneously in synchronism with 
the multiplexed time reference signals. Sampling was 
initiated at the start of an even 10 sec. interval (UT2) 
by a start pulse derived from the recorded time code. 

The sampling programs were improved between the Apollo 15 
and l6 experiments to obtain higher playback rates for the 
13 cm data (see Block Diagram I, Notes, 4d) . With this 
exception, the same computer programs were used for data 
reduction from the three sets of observations. A small 
overlap was provided between subsequent sample data tapes. 
These overlaps were carried through the remainder of the 
data processing. Tape recorder playback levels were 
adjusted for equality between left circular and right 
circular polarization using the controlled reference 
signal levels for calibration. The playback recorder 
electronics were equalized for the particular tape source 
(FRlIj.00 A, FRlliOO B, HP 3955 B) prior to sampling data from 
that source. Output from the data sampling process was 

stored on magnetic tape. The quantization level was 8 bits. 

t h A 

The i data sample generally will be denoted : data 

from the left circularly polarized antenna d. , data from 

the right circularly polarized antenna 
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Subtask 3 


Conversion to Frequency Domain 

All data have been rendered in the frequency domain 
in the form of modified complex Fourier coefficients 
( Blackman and Tukey , 1958), A sequence of data samples 
was multiplicatively weighted with a sine-squared (Hanning) 
data window and then Fourier analyzed using fast Fourier 
transform techniques. Analytically, the data samples 
were grouped, separately for each polarization, according 
to 


where i = nN 4 - j, j^N; n,N,j positive integers or zero. 
The modified complex Fourier coefficients are 


N-l 

f k = y sin2 (1 2 j + £) d j e i(2ir/N)jis (2) 

j = ° 

where i s= /-l, 0 < k <, N-l, In the Apollo data reduction 

programs, N was either 102^ or 2 0k8 (see Table Xl), 

Each set of coefficients corresponds to a time interval 

T = N*(sampling interval) (see Table II for values). 

The outputs of the transformation are the f n above, 

k 
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Subtask k 


Generation of Coherency Matrix 

The coherency matrix ( Born and Wolf , 1959) was determined 
directly from the f£ by forming the summed products 


L 


I\sT 


l 

n=l 


f n 
1 k 


f n* 
2 k 


x 

nal 


f n * 
lk 


f 

2 k 


L 



n=l 


( 3 ) 


where * denotes complex conjugate. Each sum is a function 
of the frequency index k. In order to compensate for the 
non-uniform effects of the receiving system filters, we 
further form 


M 


q 2 

l q k 


- < l4* 


ni2 


/_ 
n=l 


I q k' 


(4) 


a2 

k 


M 


K 2 q 


- I I aI 


n=l 


n n *2. 

where q = f for periods with signal absent. The q 
k k K 

are the receiving system output power spectra for a noise 

input. The upper bound M was chosen to reduce the fluctua- 
? ~ Ij. 

tions in q to a small value. Typically M < 10 , for 

p k 

which Sr is determined to approximately one percent. The 
k 

coherency matrix, corrected for receiver power transfer 
characteristics is 
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Subtask 4 (cont.) 


where 


J n "2 

l q k n=l 


I I XI 


22 .2 

2 q k n=l 


I 1 XI 


12 


Y f n f n * 

a. L ik 2 k 


l q k 2 q k 


n=l 


(6) 


J 


21 



The fractional polarization of the received signal may be 
computed directly from the J as 

™K 

4 Det J k ) 1/2 

Trace 2 J, J . 

-k 7 




Again y^_ is a function of the frequency index k. The 
output of Subtask k consists of the spectra . v . The 
parameters L and N used in the reduction of the several sets 
of observations are included in Table II, Completion of 
the J (Subtasks 2, 3* and 4) represented the greatest portion 
of the computational expense for this experiment* Additional 
results may be obtained directly from J and y * For example, 

“ K K 

the power in the polarized and unpolarized parts of the echo 
is : 


i4 




Subtask 4 (cont.) 


polarized power P^(k) “ *Y k * Trace 

( 8 ) 

unpolarized power P^(k) = (l~^ k )* Trace £ k • 


Other parameters of the echo spectra may be obtained simi- 
larly ( Born and Wolf , 1959 J Appendix 10) . 
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Subtask 5 


Correction for System Polarization Parameters 

The quantities and discussed in the previous section 
under Subtask 1| were derived directly from the sample 
data as they came from the analog tapes. As discussed 
under Subtask 4 certain corrections have been made for the 
receiver filter characteristics. However , the data were 
treated as though they were derived from perfect antennas. 
That is, the antennas were assumed to consist of a pair of 
right and left circularly polarized elements. It was fur- 
ther assumed that, with the exception of the filter correc- 
tions already applied, the gains in the two receiver 
channels were equal. In the case of the 13 centimeter 
data this assumption was very good. The isolation of the 
DSN 64 m antenna has been measured as > 2 6 db with an 
axial ratio of < 0,8 db (private communication, D. Bathker, 
JPL , 1973 ). At 116 cm the properties of the SRI 46 m 
antenna are not nearly so well known. Polarization of the 
46 m antenna was controlled principally through the con- 
straints applied to the construction of the feed system. 

The feed system consisted of a crossed dipole array of 
linear elements connected through a standard hybrid to 
obtain a circular polarization. The array elements were 
mechanically and electrically identical. The hybrid 
combiner and associated phase shift elements were adjusted 
to within 1 and 1 db of the ideal transfer function for 
such a device. Coupling between the orthogonal linear 
array elements was less than 40 db with the feed removed 
from the dish. Cross-coupling between the two circular 
polarizations, observed at the hybrid output with the feed 
in place at the focus of the dish^was l 6 db. We were unable 
to measure the axial ratio of the overall system with the 
feed in place in the dish* We estimate axial ratio of the 
116 cm system as approximately 1 db for either polarization. 
For certain received polarizations, axial ratios 




Subtask 5 (cont.) 


of this magnitude can introduce significant errors 
in the calculation of fractional polarization. 

Consequently, a correction for this uncertainty was applied 
in Subtask 5. This correction was applied to the ll6 cm 
Apollo Ik and 1 6 data only. The output tapes from Subtask 
5 preserve the original J computed earlier, 

— K 

The corrections were determined as follows (cf Tyler , 1970), 
Consider the signals arriving at the antenna terminals in 
terms of their right and left circularly polarized compo- 
nents, which we will denote E^, , respectively. The 

relationship between the arriving signals and the signals 
at the antenna terminals may be expressed as a matrix 
multiplication 


E 


E 


f 

r 

t 

1 





(9) 


where the c f s are complex and arbitrary. The matrix 
elements may be thought of as the transmission coefficients 
of the four port network consisting of pairs of antenna 
elements and terminals. Physically, the c f s may represent 
attenuation, gain, and cross-coupling. In an ideal system 

C u - °12 = C 21 - °> °22 - 1 ' 
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Subtask 5 (cont.) 


The effect of such a transformation on the coherency matrix 
J is easily shown to be 


( 10 ) 


where J is the coherency matrix of the wave associated 
with J observed at the antenna terminals. The k subscripts 
have been suppressed for convenience. However , it is 
assumed that the c*s are independent of frequency over 

t 

the spectrum of interest. Given J and the c f s, the 

original J may be recovered through an inverse matrix 
” ? 
manipulation. In the present case J is observed^ but 

the c T s are unknown. 

The c*s may be estimated from an observation of an unpola- 
rized signal ( Tyler , 1970). System noise inputs to the 
ll6 cm receiver system were used to estimate the c*s and 
perform a correction. The method was based on an experi- 
menter selection of those portions of the receiver output 
spectrum that contained only receiver noise. If it is 
assumed that the noise input is unpolarized^ then 


11 

t 

r i2 

t 

r 2l 

t 

*22 


' 11 ' 


c M c 


11 12 11 12 1 12 


°11 C 21 C 11 C 22 °12 C 21 C 22 C 12 


c *c , c c * 


c c , c c 


11 21 21 12 22 11 22 12 


c c 

1 21 22 21 


C 22 C 21 ^22^ 


11 


12 


21 


22 



(11) 
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Subtask 5 (cont.) 


where 


J 

o 

k 

T 

sys 


k T /2 
sys 


1.38 x 10 ^ (Joules/deg. Kelvin) 
system temperature. 


(12) 


The signal will appear unpolarized if 



= 0 and J 


11 



or 


c 


11 C 21 “ 


c c 

22 12 


and 



(13) 


from which the required inverse transformation can be 
obtained. In practice^ a numerical estimate of the 
correction matrix was obtained as described above. This 
estimate was then used as a starting point in a search to 

find the c T s which minimized the apparent polarization of 

? 

the corrected J^in the noise portions of the spectrum. 


The corrections in the form just described were then used 

9 

to compute the corrected J for the entire spectrum, A 

k t 1 

corrected fractional polarization y based on J was then 

K iC 

obtained. 


The output tapes from Subtask 5 contain the original J, 

t “ k 

and the new, corrected y^. A data user may easily 

recompute the original y from the J which have been pre- 

served and which are available on the tapes supplied 

(see Subtask 6, Appendix IV) . The correction factors 

? 

employed in the generation of the y are given in 

K 

Appendix VIII. 
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Subtask 5 (cont , ) 


In summary, steps in this subtask are: 

a) Read source tapes from Subtask 4 

b) Determine elements of the correction matrix 
based on minimization of the polarized part 
of the noise. 

c) Compute y 

K 

d) Generate new^tape containing the original 

J n and new y 1 „ — — 

k ’ k 

This process was applied to the Apollo ll+ and 1 6 ll6 cm 
data: 13 cm data are uncorrected, i.e., for 13 cm data 

t t 

Y, = y 1 • In subsequent steps the y were used in all 
K K k 

computations of the polarized and unpolarized parts of 
the ll6 cm echo spectrum. 
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Subtask 6 


Merging Observations with Trajectory 

t 

The output of Subtask 5, the y^ tapes^ represents the 
experimenter 1 s best estimate of the received spectra, 
averaged over the time intervals previously defined: that 
is, a complete second order description of the received 
echo signal. In the present step this data, which was 
previously processed without regard to lunar coordinates 
or other geophysical considerations, was combined with the 
Apollo command-service-module ephemeris. The emphemeris 
was first interpolated to the mid-point of the averaging 
period used in the computation of the then certain 

ancillary quantities were computed. The interpolated 
ephemeris and the derived quantities were then merged 
with the experimental data to form a basic set of 
source tapes designated JM Doptrack . After further pro- 
cessing in Subtask 10 , these tapes became the primary 
source tapes supplied with this report to the NSSDC. 


The ephemeris based quantities added to the observational 
data were : 

a) Time corresponding to mid-point of averaging 
period 

b) Predicted difference between reflected and direct 
doppler shifts 

c) Predicted echo bandwidth for a moon with rms 
slope of 0.1 

d) Angle of incidence on mean spherical moon 

e) Spacecraft altitude above mean spherical surface 

f) Speed of the spacecraft 

g) Bistatic-radar cross-section of a smooth, 
perfectly conducting, spherical moon for the 
current sp acecraf t-moon-earth geometry 

h) Normalized signal strength for a conducting moon 
and instantaneous geometry 

i) Spacecraft position in selenographic coordinates 

j) Specular point position on a mean spherical lunar 
surface in selenographic coordinates 

^3 ^ 



Subtask 6 (cont.) 


k) Selenographic latitude and longitude of space- 
craft position 

l) Doppler shift due to rotation of the earth 

m) Total doppler shift of the reflected signal 

n) Selenographic latitude and longitude of specular 
point on the mean spherical lunar surface 

o) Speed of the specular point on the mean lunar 
surface 

p) Look angles to earth in spacecraft coordinates 

q) Euler angles of spacecraft attitude in local 
horizon system 

r) Selenographic unit velocity vector of spacecraft 

s) Selenographic unit vector location of earth. 

These tapes are organized in groups of six data records , 
referred to as a data frame, corresponding to each time 
interval * A complete description of the JM Doptrack 
tape formats and contents is given in the appendixes 
(see Appendix IV)«, 
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Subtask 7 


Computation of Polarized and Unpolarized Parts, Carrier 
Suppression 

JM Doptrack tapes contain spectrally analyzed 13 cm and 
ll6 cm receiver outputs. No provision was made within the 
receivers or Subtasks 2* 3 t 5> or 6 tor removal of the 
directly propagating telemetry carrier from the echo data. 

In terms of the polarization parameters, this signal cannot 
be removed completely. However, much of the data analysis 
is based only on the low order moments of the polarized 
part of the echo. Subtask 7 computed polarized and unpo- 
larized spectra from relation (8), Subtask k, then used 
an empirically derived algorithm to remove the carrier 
signal from the polarized spectra; these data together 
with the ephemeris data described in Subtask 6 constitute 
the intermediate data set generated by Subtask 7- The 
algorithm for carrier suppression is described in Appendix V # 
Carrier suppression was also applied to unpolarized 
data. However, the carrier was largely suppressed in 
those data by coherency matrix processing, since it was a 
polarized signal. The output of this subtask, the P tapes, 
contain the experimenter * s best estimate of the power 
spectra of the polarized and unpolarized components of the 
echo signal. 
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Subtask 8 


Computation of Total Received Power, Echo Moments 

Inputs to this subtask were the polarized and unpolarized 
power spectra derived from the JM Doptrack source 
tapes (Subtasks 6, 7 ), The purpose of Subtask 8 was to 
derive numerical measures of the echo spectra. Certain 
measures of the echo spectra were corrected for predictable 
trajectory effects and converted into scientific units. 

In all cases, the designation of the echo signal location 
in a spectrum was made by the experimenter * s visually 
scanning plots of the polarized and unpolarized power spec- 
tra* The values of the k indices bounding the echo were 
input to a computer program that did the actual data reduc- 
tion. An average noise level, determined from a region of 
the spectra not containing echo, was also input to the 
computer program. The data were monitored at approxi- 
mately 30 sec. intervals and the echo limits reset to 
account for motion of the echo in the receiver passband. 

It was also necessary to occasionally reset the noise 
level as it also varied during the experiment, principally 
due to the scan of the receiving antenna across the lunar 
terminator. The derived quantities are given below: 

a) Polarized echo power — the integral of the 
polarized power spectra between the frequency 
limits set by the experimenter and above the 
system noise level. This quantity is the best 
measure of the polarized echo power received. 

b) Normalized polarized echo power — the quantity 
described in a) above, divided by the polarized 
system noise level. 

c) Unpolarized echo power — the integral of the 
unpolarized power spectra between the frequency 
limits set by the experimenter and above the 
system noise level. This quantity is the best 
measure of the unpolarized echo power received. 

d) Normalized unpolarized echo power — the quan- 
tity described in c) above, divided by the unpola- 
rized system noise level. 


2k 



e) Equivalent area bandwidth — the bandwidth of 
the polarized echo signal between the frequency 
limits set by the experimenter and above the 
polarized system noise level as determined by 
the ratio of the total polarized echo power to 
the peak polarized echo power. 

f) Normalized absolute moment bandwidth — the band- 
width of the polarized power echo as described in 
e) above, computed from a gaussian equivalent 
absolute moment, and divided by e). 

g) Normalized second moment bandwidth — the 
bandwidth of the polarized power echo as described 
in e) above, computed from a gaussian equivalent 
second moment, and divided by e). 

h) Centroid of echo spectrum — the centroid of the 
polarized power echo as described in e) above. 

i) RMS slope — the rms slope of the lunar surface 
inferred from the value of e) above and the pre- 
dicted bandwidth for an rms surface slope of Q.l , 
using linear interpolation. 

Formulas for the computation of the above quantities and a 
brief explanation of their use may be found in Appendix VI. 
The trajectory data added in Subtask 6 were retained through- 
out this subtask. A simple correspondence between the out- 
put of this subtask and the JM Doptrack tapes was maintained 
through inclusion of ephemeris data in both data sets. 



Subtask 9 


Final Editing of Data 

Final editing of the output from Subtask 8 resulted in a 
set of reduced data records designated Integral tapes. 

These tapes constitute the second form of data supplied 
to the NSSDC. 

Final editing of data included the following steps: 

a) Addition of hand-scaled bandwidths as a partial 
independent check on Subtask 8 (see Subtask 10 ) . 

b) Addition of spacecraft antenna gain in the direction 
of specular reflection. 

c) Notation of operational or data processing changes 

i) Polarized band bad 

ii) Unpolarized data bad 

iii) Change in polarized noise level 

iv) Change in polarized k indices for echo limits 

v) Change in unpolarized noise level 

vi) Change in unpolarized k indices for echo 
limits 

vii) Change in system gain 

Caution: Some integral data, primarily at the beginning 

or end of a transmission, have been deleted. Large over- 
laps in the 13 cm data caused by the use of two analog 
tape recorders have also been deleted. In all other cases 
bad data, for example when interference is present, are 
flagged on the edited Integral tapes as described in c) 
above (see Appendix VI, Ei+T ) • The data contained on the 
Integral tapes are still the experimenter* s best estimates 
of the values. However, in the case of flagged data, that 
estimate may be very poor. No flagged data should be used 
without examining the spectra on the corresponding JM Dop- 
track tapes. For example, in the ll6 cm data occasional 
interference produced marked increases in the apparent 
polarized echo power, but evidently left the unpolarized 




Subtask 9 (cont.) 


power unaffected, A flag for polarized power thus may 
also cast suspicion on the unpolarized power. The experi- 
menters have evaluated these cases and indicated their 
opinions accordingly in the data. Other individuals may 
arrive at different conclusions. It is primarily for 
this reason that the JM Doptrack/Integral tape 
f rame-to-record correspondence (through ephemeris data) 
has been maintained: reduced data records on the Integral 
tapes may, if questioned, be re-evaluated from the JM 
Doptrack source tapes, A complete description of the 
Integral tape formats and the flags is given in Appendix VI. 



Subtask 10 


Displaying, Scaling, and Copying the Data 

Output from Subtask 6, the JM Doptrack tapes, is used in 
miscellaneous programs in preparation for data analysis 
and dis ti'ibution. An important step in producing the 
final Integral tapes (Subtask 9 ) took place here* Polari- 
zed power spectra obtained from JM Doptrack tapes are com- 
puter plotted and visually examined to determine specifically 

a) proper receiver operation 

b) proper receiver tuning 

c) presence of interference. 

Such plots are also used to determine the one-half power, 
hand-scaled bandwidths added to the data set in Subtask 9* 
This bandwidth is determined by measuring the width of 
the polarized echo spectrum at a point one-half the 
distance from the apparent system noise level to the mean 
echo peak* The measured distance is scaled by the 
appropriate factor to determine the width in Hertz. 

Such measures can be related to the rms slope of the lunar 
surface from the quasi-specular scattering theory 
(see Appendix Vi). The hand-scaled values are used to 
verify the machine algorithms used in computing lunar 
rms slopes and as a simple means of quickly estimating 
the slope. 


Finally, this subtask included copying of the JM Doptrack 
tapes for shipment to the NSSDC* These tapes were not 
edited. 
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Subtask 11 


Utility Routines 

Integral tapes have been organized as a sequence of records 
describing the data and experimental geometry at successive 
instances of time. It is then a simple matter to determine 
any set of variables from this tape with time as a parameter. 
We suggest that data users consider this data as a set of 
dependent functions parameterized in time and hope that such 
a presentation is found useful. 
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Stanford Apollo ll6 cm 


(VHF) Bist atic~Radar Receiving System: 
Major Elements 
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0 Hz CHANNEL I 5000 Hz 


Apollo Bistatic-Radar Receiving System 
ll6 cm Uncorrected Frequency Response 
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1 I 

0 Hz CHANNEL 1 21500 Hz 


Apollo Bistatic-Radar Receiving System 
13 cm Uncorrected Frequency Response 
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Appendix IV 

JM Doptrack Tape Formats 


This appendix describes the tape formats for the JM Doptrack Tapes 

generated in Subtask 6 of the text. Tapes are 9 track, binary in XDS 

<y 

Sigma 5 machine images. These tapes contain the output of the polari- 
meter in Subtask 4 , the corrected fractional polarization obtained from 
Subtask 5, MSC trajectory data, and certain ancillary quantities compu- 
ted at Stanford. All records within the given file are the same length, 
There may be more than one tape per file. Files are identified by a 
header record which contains a brief description of the file contents. 


This appendix 

describes 

the tape organization, the file organization, 

and the record 

1 formats for the JM Doptrack 

tapes 

. Definitions of the 

tape contents 

are either 

given here or 

described 

by reference. 

A. File Organization 





File 

No. Tapes 

Contents 


Record Length 

1 

1 

Apol lo-l 4 

116 

cm 

514 words 

2 

1 

Apollo-lU 

13 

cm 

514 words 

3 

1 

Apollo-15 

116 

cm 

1026 words 

h 

2 

Apollo-15 

13 

cm 

514 words 

5 

1 

Apollo-l6 

116 

cm 

1026 words 

6 

2 

Apollo-l6 

on 

r-t 

cm 

514 words 


B. File Organization 

1 . Header Record 

2 . Data Record 1 

Data Record 2 

Data Record 3 

Data Record 4 

Data Record 5 

Data Record 6 


Observational Data 


Ephemeris Data 

* 

* 

< Many Data Frames > 


Data Frame 


3 . EOF (End of File Mark) 


a 


See Appendix XI for a description of Sigma 5 machine images; 



Appendix IV (cont.) 


If more than one tape was needed to contain the JM Doptrack 
file, the file was continued across the end of the tape to the 
next tape without interruption; there is only one header record 
per file, at the beginning of the first tape of the file* 

Once again, even though the files are physically nothing more 
than a continuous string of records, all of equal length, it 
is helpful to think of them (see figure below) as being composed 
of a header record followed by groups (frames) of six records* 
The ephemeris data record concludes the frame and, since each 
JM data frame is a short-time average, the ephemeris data have 
been calculated using the midpoint of the frame as the instan- 
taneous time reference point* 


JM Doptrack File Containing N Data Frames 


Rec. No. : 

1 

2 

3 

k 

5 

6 

7 

Header 

Record 

Data 
Rec, 1 

Data 
Rec , 2 

Data 

Rec - 3 

Data 
Rec . 4 

1 

ES^BB 

Eph. 
Rec. 6 

Data 
Rec. 1 


1 





l 

1 




1 

DATA FRAME 1 




8 

9 . 

10 

11 

12 

13 

Ik 

15 

Data 
Rec. 2 

Data 
Rec. 3 

Data 
Rec , b 

Data 
Rec , 5 

Eph* 
Rec * 6 

Data 
Rec. 1 

Data 
Rec , 2 

Data 
Rec * 3 

1 1 — . — 


DATA FRAME 2 DATA FRAME 3 


6n-i 6n 


• 

9 

9 

Data 

Eph. 

EOF 




[ Rec. 5 

Rec . 6 


DATA FRAME N 


















Appendix IV (cont.) 


C. Header Record Formats 


D. 


Word No . 

1-42 

43 

44 
45-46 

47-48 

49 

50 


Contents 

Alphanumeric file identifier 

Day of year on which data were 
collected (January 1 = l) 

Year data were taken (Gregorian) 

Julian Ephemeris Day at 00:00 UT2 
on the day the data were taken 

Julian Ephemeris Day of reference 
epoch 


of data averaging frame 


Number of data records following 
this header record (Number of 
data frames = number of data 
records divided by 6) 

51-514/1026 Not used 

Data Frame Formats 

1# Data Record Organization 


Record No , 
1 
2 

3 

4 

5 

6 


Contents 

J n< k ) 

J 22< k > 

Real part of J.^k) 
Imaginary part of J^fk) 
Y (k) 

Spacecraft ephemeris and 
ancillary data 


I 

R 

DPR 


Alphanumeric 

Integer 

Real 

Double Precision Real 


Units 

Of 

Machine Type 

— 

A 

(days) 

I 

(year) 

I 

(days) 

DPR 

(days) 

DPR 

(sec) 

R 

— 

I 


Note: 

Il6 cm 

•V k > 

data 

- Left circu- 
lar polariza- 
tion 

J 22< k > 

- Right circu- 
lar polariza- 
tion 

13 cm data 

J-^(k) - Right circu- 
lar poliriza- 
tion 


J^fk) - Left circu- 
^ lar polari- 

zation 
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D. 


2 . 


Format Records 
Record No . 

1-513/1025 

514/1026 


1-5 (all machine type real ) 05 
Contents 

Data described under D.l 
above 

Meaningless 


3- Format Record 6 (all machine type real ) 


Word No. 

Contents 

Units 

1 

Meaningless 


2 

UT2 at midpoint of frame 

(sec) 

3 

Reflected doppler minus 
direct doppler 

(Hz) 

4 

Predicted bandwidth for 
rms surface slope of 0.1 

(Hz) 


5 

6 


Angle of incidence (deg) 

Spacecraft altitude (mean 
lunar radius assumed 
1736 km) (km) 


7 Spacecraft speed (m/sec) 

8 Radar cross section pre- 

dicted for smooth conduc- 
ting moon (dim) 


9 

(Radar cross section)/ 

0 


(received power) 

(m yw) 

10 


Components of seleno- 


11 

Y 

graphic unit position 
vector of spacecraft 


12 

Z 

J 

location 

(dim) 

13 

x' 

Components of seleno- 


i4 

Y 

graphic unit position 
' vector of specular 


15 

Z 

point location 

(dim) 

16 

Selenographic latitude of 
subspacecraft position 

(deg) 

17 

Selenographic longitude of 
subspacecraft position 

(deg) 

18 

Component of doppler shift 
due to earth rotation 

(Hz) 


°^)ue to a hardware problem in the data processing, words 1-15 are zero 
in the Apollo-l 4 , ll6 and 13 cm, and Apollo-15 13 cm files. 
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D. 3* Format Record 6 (all machine type real ) (cont.) 


Word No. 


Contents 

Units 

19 

Total doppler shift of 
reflected signal 

(Hz) 

20 

Selenographic latitude of 
specular point 

(deg) 

21 

Selenographic longitude of 
specular point 

(deg) 

22 

Speed of specular point 



on 

the lunar surface 

(M/sec) 

23 






y Vehicle look angles 


2k 

Pe 

to earth 

> 


25 

0 ' 



26 

7 

Euler angles of local 
K horizon coordinates 

(deg) 

27 




28 

X 

Spacecraft selenographic 


29 

Y 

y unit velocity vector 

(dim) 

30 

z. 



31 

X 1 

Selenographic unit 


32 


vector from center of 


Y 

Y | 

moon to center of 


33 

A 

earth 

(dim) 


34“5l4/1026 Not used 


Note? Data do not always progress uniformly in time. Occasionally, 

data frames will reverse in time for one frame, and then continue 
forward. This effect is caused by the sampling procedure in 
which deliberate overlap was inserted. Time on data is correct. 


E ♦ Definition of Contents 

1. Data — defined under Subtasks 4 and 5 of the text. 

2 * Ephemeris and Ancillary Data — described in greater detail in 
Appendix VI. 



Appendix V 

Carrier Suppression Algorithm 

This appendix describes the algorithm used to remove the direct 
signal from the polarized and unpolarized parts of the bistatic-radar 
echo spectrum (see Subtask 7 ). This algorithm was used on each frame 
of polarized and unpolarized data which the JM Doptrack tapes yield 
according to relation (8), Subtask k* Two slightly different procedures 
were used for the ll6 cm and 13 cm data. These differences constituted 
the only difference in procedure between the lib cm and 13 cm data. 

Both data sets were processed by the same computer program, with differ- 
ent program branches for the two cores. The direct signal observed at 
13 cm was free of spurious sidebands to the level of our observations. 

At llo cm, the direct signal contained two weak sidebands symmetrically 
displaced approximately 20 db below the direct signal. Consequently, 
during periods of strong direct signal at the ll6 cm wavelength, it was 
also necessary to correct for the presence of these sidebands. The pro- 
cedures described below were developed empirically but were found 
to give good results. The approach was to find the maximum of the power 
spectrum and to assume that this maximum represented the direct signal. 
This assumption was tested by determining the height of the maximum with 
respect to the fluctuations in the spectrum in the immediate vicinity 
of the maximum. If the maximum exceeded the fluctuation criteria then 
an interpolation procedure was used to provide a smoothed estimate of 
the spectrum. In the case of the ll6 cm data, sideband suppression was 
achieved by reducing the local maximum on either side of the direct 
signal by an amount proportional to the strength of the direct signal. 

A detailed summary of this procedure is given below. 


Initial Test 


1, Find the absolute maximum of all the data, MAX (see figure 
which follows). 


2. Find the adjacent local maxima, MAX12 and MAX34. 

3. Test MAX to determine whether or not it represents the direct 
signal . 


i) 


SUM = MINI ± MINg + M AX 12 

4 


+ 


MIN3 4- min4 


+ MAX 34 



8 


39 


8 
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Initial Test (cont.) 

ii) GRIT = 1* 2-SUM 

iii) ( MAX ♦LT. GRIT) = 0 > no direct signal present , terminate proce- 
dure; go to next frame, initial test 1 # 

(MAX. GE. CRIT) =£> direct signal present, continue. 

4* Check data type. If data is ll 6 cm then go to step 8 , otherwise 
continue with step 5 . 


MAX 



13 cm Procedure 

5 . Test for a smooth direct signal. 

IF (2*0*MIN1.GE*MIN2) 

and =£> smooth direct signal* 

( 2 . 0 *min 4 .ge .min 3 ) ' ' “ 

If the direct signal is smooth then all values between MIN2 and 
and MIN3 which are greater than CRITC are set equal to CRITC 
and the procedure is terminated; initial test 1 is then begun 
on the next frame. If the test for smoothness is failed, 
then continue with step 6 , 


6 . This procedure is applied only if the direct signal is spread 
through several frequency analysis bins. Determine the first 
minimum on either side of the direct signal which 
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Initial Test (cont.) 

satisfies the following condition: 

(MIN(l) ♦LT.MIN(l-l)* 2*0) Where I indexes the left side of 

the direct signal 

(MIN( J) *LT *MIN( J+l) # 2* 0) and J indexes the right side of 

the direct signal (see fig, below) 

If the conditions are not satisfied in 15 minima^ choose the 
15 th minima* Denote the extreme minima I, J, MIN(l) and 
MIN(j) locate the extent of the spread direct signal. 

Denote the minima adjacent to the direct signal I, J. Then 
replace the data D( ) with D f ( ) calculated as follows: 

D«(L) = D(L) - -(*) - D (j . ) • (L-l) ; I + 1SLSI-1 

1 - 1 

D' (L) = D(L) - ~ D ( J ) , (L-J) ; J + 1 £ L <; J - 1 

J - J 

D»(L) = D(I) + Si ?) - D ( J ) ( L -i); 

J - I 

7. Terminate procedure; go to next frame, initial test 1. 
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Il 6 cm Procedure 

8 . Find the largest local maximum within 60 Hz of the direct 
signal on first the left side and then on the right side 
of the direct signal. Designate these maxima MAXL, MAXR, 
respectively. 

9 . Determine relative magnitude of the direct signal. 

CRITSB = 200.0' SUM 

If (MAX. GE .CRITSB) then the direct signal is sufficiently 
large to require compensation for the sidebands. Replace 
MAXL and MAXR by the average value of their respective adja- 
cent minima. 

10. If (MAX.LT. CRITSB) the direct signal is small. Subtract a 
constant from the maxima MAXL and MAXR, defined above. The 
constant is 

ESUB = (MAX - SUM)/500.0 

The subtraction is performed only if the data value is greater 
than ESUB. 

11. Replace MAX by the average of MIN2 and MIN3. 

12. Terminate procedure and go to next frame, initial test 1. 



Appendix VI 
Integral Tape Format 

The Integral tape contains reduced data records generated in Subtask 
1 » 9 9 (see Block Diagram I ) . This appendix gives the detailed for- 

matting of that tape and describes all ancillary computations. 

A. Tape Organization 

The table below gives the file contents and data record length for 
the Integral tape. 


File No . 
1 

2 

3 

4 

5 

6 


Contents 

Apollo-l 4 ll6 cm 
Apollo-l 4 13 cm 
Apollo-15 ll6 cm 
Apollo-15 13 cm 
Apollo-l6 ll6 cm 
Apollo-l 6 13 cm 


Record Length 
50 words 
50 words 
50 words 
50 words 
50 words' 

50 words 


The tape is 9 track, 800 BPI, binary in XDS Sigma 5 machine images. 
See Appendix XI for a description of Sigma 5 machine images. 


B. File Organization 

1 . Header Record 

2 . Data Record 


<many data records^ 


3 . (EOF) End of File 

4 . Header Record 


5 * - 
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C. Header Record Format 


Word No* 


Contents 


Units 


Machine Type 


a 


1-42 

43 

44 
45-46 

47-48 

49 

50 


Alphanumeric file identifier — A 

Day of year on which data were 
collected (January 1 = l) (days) I 

Year data were taken (year) I 

Julian Ephemeris Day at 00:00 UT2 
on the day the data were taken (days) DPR 

Julian Ephemeris Day of reference 
epoch (days) DPR 

Time increment between center- 
point of data averaging frame (sec) R 

Number of data records follow- 
ing this header record — I 


D„ Data Record Format (all machine type real ) 


Word No. 



Contents 

Units 

i 



Meaningless 


2 

UT2 

at midpoint of frame 

( sec) 

3 

X 1 


Components of seleno- 


4 

Y 

► 

graphic unit position 
vector of spacecraft 

(dim) 

5 

z 


location 


6 

x 1 


Components of seleno- 


7 

y 

, 

graphic velocity unit 
vector 

(dim) 

8 

. Z J 




9 

Speed: Magnitude of space- 

craft velocity vector 

(m/sec 

10 

X 1 


Components of seleno- 


n 

Y 

> 

graphic unit vector from 
center of the moon to 

(dim) 

12 

Z 


center of the earth 



( y 

A - Alphanumeric 
I - Integer 
R - Real 

DPR - Double Precision Real 
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D. 


Data Record Format (all machine type real) (cont.) 


Word No . 

13 

14 

15 


Contents 


X 

Y ► 
Z 


Components of seleno- 
graphic unit position 
vector of specular 
point location 


Units 


(dim) 


16 

17 

18 


e 

♦ k 

0. 


Euler angles of space- 
craft attitude and 
local horizon frame 


(deg) 


*»• J " [(*)] [(*)] [(*)] [ a i ( 


cw rotation look- 
ing in + axis 
direction 


19 

20 

21 


22 

23 

2k 

25 

26 

27 

28 

29 

30 

31 

32 

33 


a g Vehicle look angles to specular 
f point (deg) 

p s 

6 Angle between plane of inci- 
dence and plane containing 
both the vehicle x axis and 
direction vector to specular 
point 


Q- ) 

e t Vehicle look angles to earth (deg) 
P eJ 

Selenographic latitude of space- 
craft position (deg) 


Selenographic longitude of space- 
craft position 


(deg) 


Selenographic latitude of specu- 


lar point (deg) 

Selenographic longitude of 

specular point (deg) 

Angle of incidence (deg) 

Instantaneous speed of specu- 
lar point on lunar surface (m/sec) 

Predicted bandwidth for rms 
surface slope of 0,1 (Hz) 

Reflected doppler minus direct 
doppler (Hz) 

Total doppler shift of reflec- 
ted signal (Hz) 

Component of doppler shift 
due to earth rotation 
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D* 


Data Record Format (all machine type real ) (cont.) 


Word No. 

Contents 

Units 

34 

Altitude of spacecraft above 
lunar surface: Radius of the 

moon assumed to be 173 6 Km 

(Km) 

35 

Radar cross-section predicted 
for smooth conducting moon 

(dim) 

36 

(Radar cross-section)/ 
(received power) 

(m 2 /w) 

37 

Polarized power 

(arb) 

38 

Normalized polarized power 

(° k ) 

39 

Unpolarized power 

(arb) 

40 

Normalized unpolarized power 

(° k ) 

41 

Equivalent area bandwidth 

(Hz) 

42 

Normalized absolute moment 
bandwidth 

(dim) 

43 

Normalized second moment 
bandwidth 

(dim) 

44 

Centroid of the echo spectrum 

(Hz) 

45 

RMS slope inferred from equiva- 
lent area bandwidth 

(deg) 

46 

Spare if value equals zero, 
otherwise handscaled one-half 
power echo bandwidth 

(Hz) 

47 

Data validity flag 


48 

Spare if value = 0* other- 
wise value of spacecraft 
antenna gain in a , p , 
direction (see word no* 19 , 

20 ) 

(dim) 

49 

Not used 


50 

Data record sequence number 



Note: Data do not always progress uniformly in time. Occasionally , 

data records will reverse in time for one record, and then 
continue forward. This effect is caused by the sampling 
procedure in which a small deliberate overlap was inserted* 

Time tags on data are correct* Overlapping data correspond 
to the same time interval but different sampling passes. 

Slight differences arise from variation in exact times averaged* 


k6 
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E. Data Parameter Definitions 

The remainder of this section defines the contents of the data 
records described just above. The individual subsection numbers corre- 
spond to the word numbers in Section D (Data Record Format). If a 
particular quantity has been described at length elsewhere, a reference 
will be given. Otherwise, the quantity is defined here. 

1. Meaningless 

Time is the OT 2 in seconds at which the data were taken. This time 
corresponds to the mid-point of the averaging interval, as described 
under Subtask 4 * All trajectory parameters have been interpolated 
to this time, so that geometrical quantities correspond to the loca- 
tion of the specular point on the mean lunar surface at the middle 
of the averaging interval. 

5 # S elenographic Unit Position Vector is defined with respect 
to the lunar surface. This vector is a unit vector directed from 
the center of the moon for the instantaneous location of the space- 
craft. The X, Y, Z directions are defined as follows: 

X s s mean earth direction 
Y = mean direction of the following limb 
Z =r north polar direction. 

The sel enographic coordinates were obtained by rotation from the 
selenocentric geo-equatorial units of the epoch given in the 
header record. Procedures are described elsewhere (Tyler, 1968 ). 

6 , 7, 8 - Selenographic Unit Velocity Vector is a unit vector in the 

direction of the spacecraft velocity. The coordinate system is the 
same as that given in items 4 , 5 above. 

9. Speed is the magnitude of the spacecraft velocity vector. 

10, 11, 12* Selenographic Unit Vector to Earth is a unit vector giving 
the direction from the lunar center of mass to the center of mass to 
the earth in the selenographic coordinate system described in the 

3 , 4, 5 above. 
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13 , l 4 , 15. Selenographic Unit Position Vector of Specular Point 

is a unit vector from the center of mass of the moon to the loca- 
tion of the specular point on the mean spherical lunar surface. 

For this computation the lunar radius was taken as 1736 km. The 
specular point is the location on the mean lunar surface where the 
angles of incidence and reflection are equal (cf Tyler, 1968). 

l6, 17, 18 . Euler Angles of Spacecraft Attitude connect the space- 
craft altitude with a local horizon reference frame. Both the 
Local Horizon system and the Euler angles are defined in Appendix IX. 

19 > 20 . Look Angles to Specular Point are in vehicle polar coordinates. 
These quantities are defined in Appendix IX. 

21 . Plane of the Vehicle , the orientation with respect to the plane of 
incidence, is given by the angle 6. This quantity is necessary to 
define the vehicle attitude with respect to the plane of incidence. 
The angle 6 is defined in Appendix IX. 

22 , 23 * Look Angles to Earth are the vehicle polar coordinates of a 
unit vector in the earth center of mass direction. These quantities 
were computed in the same manner as items 19, 20 above using the unit 
vector to earth. 

2 k. Selenographic Latitude of Spacecraft Position is the selengraphic 
latitude of the subspacecraft point computed from the Z component 
of the unit vector given in items 3, k, and 5. 

25 - Selenographic Longitude of Spacecraft Position is the selenographic 
longitude of the subspacecraft position computed from items 3, 1+, 
and 5 according to astrometric convention, western limb of the moon 
leading. 

2 6. Selenographic Latitude of Specular Point is the selenographic 

latitude of the specular point on a mean spherical lunar surface 
computed from item 15. 

27 * gelenographic Longitude of Specular Point i s the selenographic 

longitude of the specular point on the mean spherical lunar surface 
computed from items 13, l 4 , 15. 
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of Incidence is the angle of incidence on mean spherical lunar 
surface at the specular point (cf Tyler, 1968). 

29 - Speed of the Specular Point is the speed with which the instan- 
taneous specular point moves across the mean lunar surface (Tyler, 
1968). 

30 * Predicted Bandwidth is the one— half power spectral width predicted 
for an rms surface slope of 0.1, based on the instantaneous angle 
of incidence and specular point velocity. Computation is after Fjeldbo 
(1964)* also described in Tyler , (1968). Fjeldbo gives a theoretical 
expression for the 1/2 power echo bandwidth: 

Af = 4(21n2) 1/2 Jocos^, 

* d o 

where v = velocity of the specular point on the mean lunar surface, 
s 

X tss wavelength of the radiation (either ll6 cm or 13 cm), 4 > = angle 

h 

of incidence at the specular point, and the quantity o = the mean 
lunar rms slope* d o 

The quantity Af is the 1/2 power bandwidth predicted for a gaus- 
sian spectrum. Such a spectrum would result from a gently undulating 
surface with gaussian autocorrelation function, 

3 i* Difference Between Reflected and Direct Doppler Shifts is the pre- 
dicted frequency difference between a wave reflected from the specu- 
lar point and the signal traveling directly from the spacecraft to 
earth. Sign convention is such that the difference is positive for 
a reflected doppler shift greater than the direct doppler shift. 

32 * Doppler Shift is the total doppler shift expected from the reflected 
signal. Computation of this doppler shift included spacecraft motion 
and the earth 1 s rotation, but did not include the rate of change 
of distance between the earth and the moon. 

33 . Doppler Due to Earth 1 s Rotation is the component of the observed dop- 
pler shift due to the earth* s rotation, for a signal arriving from the 
direction of the moon. 
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34 . Altitude of the spacecraft above the lunar surface has been compu- 
ted assuming a lunar radius equal to 1736 km. The magnitude of the 
spacecraft radius vector from the lunar center of mass is obtained 
by adding the contents of word no. 34 to 1736 km. 


35 . Normalized Bistatic-Radar Cross-Section is the bistatic-radar 
cross-section of a smooth conducting sphere of the same radius 
and relative geometry as the moon. Following Fjeldbo (1964) this 
cross-section is given by 


°B = 


4jt cos < t > 


(cos , + 121) ( 


1 + 


2 d cos <(> 
or 


where 


= distance from transmitter to the center of the moon 

£ 

R = lunar radius (1.736 x 10 m) 


<t> = angle of incidence (cf item 28) 


d ” distance the transmitter to the specular 

point on the mean lunar surface. 


3 ^' (Radar Cross Section)/ (Recei ved Power) is a multiplicative .constant 
relating instantaneous geometry and received power to surface reflec- 
tivity. 


radar cross-section _ 
received power 


( 4 m:) 2 R* 

A P T °T CT B 


where 


= distance from transmitter to center of the moon 
- distance from receiving site to center of the moon 
A = effective aperture of receiving antenna 
P T = transmitted power 

G^ = transmitting antenna gain in specular point 
direction 

cTg = bistatic-radar cross-section for a perfectly 
conducting moon* 

For convenience t this expression was evaluated with the following 


numerical values for the quantities above: 
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= instantaneous value from MSC trajectory 

a* ss instantaneous value from item 35 a h° ve 
B g 

R = 1| x 10 m 

2 P 

A a 0,5 (22.5) * 

G_ = 1 
T 

P T = 2-5 w. 

These values give only order of magnitude results for this experi- 
ment , 

37 . polarized Power is the experimenter’s best estimate of the polari- 
zed component of the received echo total power. Extraction of the 
polarized power is discussed elsewhere (cf Subtask 8 ). 

Denote the polarized power spectrum P^(k) , Consider the figure 
below. Polarized power is determined from 


B, 


p = X ( p P ( k ) 


k= B, 


The P (k) is a polarized power spectrum. In the determination of 
P 

P, the signal limits and n were selected by the experimenter. 

The quantity n was chosen on the basis of 


B, 

c 

l 

b! 


( p p (k) -n) ~ 0 


where represent spectral limits containing no echo signal > 

and n represents a sequence of trials of n. The limits 

were varied as is necessary to follow the changing echo signal, 
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38. Normalized Polarized Power is the quantity contained in item 37 
divided by the average power spectral density of the system noise 
level. This quantity has been discussed in detail elsewhere 
(cf Subtask 8) , 

Using the notation introduced under item 37, 


normalized polarized power = P/n, 
where P and n have the same meaning as above. 

Note: P/n is extremely sensitive to the choice of n P Thus, polarized 

power is considered the best overall measure of received polari- 
zed echo power. But P/n provides the only method, through the 
measure of system temperature, of obtaining an absolute power 
calibration. Similarly, the value of n may be determined from 
the ratio of polarized power to P/n, so that the variations and 
system temperature and/or gain may be estimated, 

39 ' Unpolarized Power is the analogous quantity to item 37 y for the 
unpolarized power spectrum. The unpolarized power was obtained in 
a manner similar to that used to compute polarized power. Letting 
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P y (k) represent the unpolarized power spectra, and referring- to the 
figure below, the unpolarized power is given by 

*2 

U =I ( p u ( k ) “”)» 

where the tildas refer to the values of B and n used in the unpola- 
rized power spectrum. In general, the limits for the polarized and 
unpolarized echoes were different, as was the value of the 
system noise level. The difference in frequency limits arose from 
the difference in the spectral distribution of the unpolarized power; 
the difference in system noise temperature arose from the signal 
processing used to separate these quantities (cf Subtask U, 5)* 

The value of n was chosen in a manner similar to that of n in item 

37* In some cases, it was not clear that all the unpolarized power 
is contained in the receiver passband. In this event, or B Q was 
set equal to the upper or lower frequency limit as appropriate. 
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40. Normalized Unpolarized Power is the analogous quantity to item 38, 
for the unpolarized, power spectrum* The normalized unpolarized 
power is defined as 


normalized unpolarized power = U/n 

where the symbols have the same meaning as in item 39 * Comments 
given under item 38 are also germane to normalized unpolarized 
power . 

* 

4l. Equivalent Area Bandwidth : the spectrum of a bistatic-radar echo 
from a well behaved surface may be written as (Fjeldbo, 1964 ): 


s ( f ) = e COS <t> (^p)"] =e -f 2 /2o- 2 

o / f 

where 

f = frequency measured from the centroid of the echo spectrum 

v = speed of the specular point across the mean lunar surface 

X = wavelength 

<f> ^ angle of incidence , and 

h /d = unidirectional rms slope . 

00 


The three machine calculated bandwidths, i.e., the equivalent area 
bandwidth, the absolute moment bandwidth, and the second moment band- 
width, provide three quasi— independent methods of determining the 
spectral width of the received echoes. The equivalent area bandwidth 
provides a standard result that is not particularly sensitive to the 
gaussian, or non-gaussian nature of the echo spectrum* The absolute 
moment bandwidth and the second moment bandwidth so emphasize depar- 
tures from gaussian because of the increasing importance given to the 
wings of the spectrum. 

Equivalent area bandwidths have been used to determine lunar rms 
slopes. The absolute moment bandwidths and the second moment bandwidths, 
when normalized by the equivalent area bandwidth, give a sensitive 
measure of the departures of the echo spectra from the gaussian conditions. 
RMS slopes derived from these measures are termed "gaussian equivalent 
slopes” in that they would correspond to true surface conditions for a 
surface with gaussian statistics and a gaussian autocorrelation function 
with the same equivalent widths. A more complete description of the 
lunar slopes requires additional analysis (e.g., see Parker and Tyler , 1973)* 
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Solving for rms slope in terms of measured values of standard 
deviation, a, of an experimental spectrum yields 

A 

h /d « 2 . 

° ° 2(v s A ) cos <t> 

Thus, the rras slope may be readily determined from an experimental 
curve in terms of the e ^ width of that curve. The equivalent 
area bandwidth is a measure of a based on an equivalent rectangular 
spectrum of the same area as the experimental spectrum. This width 
is computed as 


i v k > 

fi 

max P (k) 



c x < k < c 2 . 


Referring to the figure which follows, the quantity a = a if the 

G3 

observed curve is gaussian and noiseless. For non-gaussian data 

ct is still a measure of the bandwidth, although the interprets- 
e a 

tion must be modified, RMS slopes determined from a and the 

ea 

expression immediately above will be referred to as equivalent 

area slopes. The quantity a is the equivalent area bandwidth* 

ea 
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^2. Normalized Absolute Moment Bandwidths are based on an equivalent 

value of a computed from the absolute moment of the data. That is, 



for a gaussian spectrum the equivalent area moments and the absolute 
moments will be equal 


A A 



where the symbols have the same meaning as in item 4l# 
normalized absolute moment bandwidth is 


a /a . 

am ea 


For a gaussian echo spectrum this ratio will be unity. 

Normalizec ^ Second M o ment Bandwidth ; the second moment bandwidth is 
also based on gaussian equivalence. This bandwidth is defined as 
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V P ( k ) (k-k ) 2 


k=C! 


sm 


I p p< k > 




For a gaussian echo spectrum 


°sm ^am ^ea 


l p „( k )- k 


k=C 


k = 


"2 

1 V k > 


k=C 


where the symbols have the same meaning as under items 42 and 4l , 
The normalized second moment bandwidth is given by 



/a 


ea. 


Again, departures of this ratio from unity are indicative of a 
non-gaussian received echo spectrum 


44# Centroid of the Echo Spectrum ; the centroid of the echo spectrum is 
defined in the standard way; 


k = 


i p P «- k 

k - C l 


I V k ) 

k=C 


where the symbols have the same meaning as in 4l , 42 and 43. The 
values of C^, are held constant throughout items 4l, 42, 43 and 
44. The echo spectrum centroid is used in the computations of the 
absolute moment bandwidths and the second moment bandwidths. It is 
also useful for estimation of the observed doppler difference 
(cf item 3l) * However, the centroid has no direct meaning in an 
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absolute sense because the transmitter frequencies from the space- 
craft are not known precisely. 

45. RMS Slopes are obtained from the equivalent area bandwidths accord- 
ing to 


h <j 

o ea 

d c " 2(v g A) cos * 


and 

unidirectional rms slope (deg.) s? tan 1 (h^/d ), 
where the results are expressed in degrees for convenience. The 

/N. 

quantities v , 4> and cj have been defined in item 1+1. 

s'' ea 

46 . Handscaled 1/2 Power Echo Bandwidths : handscaling of polarized 

echo bandwidths has been discussed previously under Subtask 10. 

A 1/2 power echo bandwidth is defined as the quantity 1f B" in the 
figure which follows* In practice, this value was obtained by 
handscaling data from plots. The handscaled bandwidths provide a 
rapid, simple technique for the evaluation of lunar rms surface 
slopes. The handscaled bandwidths were also used as a controlling 
parameter in development of automated techniques for obtaining 
echo bandwidths. RMS slopes may be obtained from the handscaled 
bandwidths by 

rms slope = 5*7° x B 

(predicted 1/2 power bandwidth for rms slope of 0.1) 

where the result is given in degrees for convenience. Predicted 1/2 
power bandwidths for 0.1 rms slopes are given in word no. 30 of the 
integral tape data records. Hand measurements were made only for 
selected segments of the data. 
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47* * the word contains a seven level binary code that indi- 

cates data quality* A bit in the "1 ,T condition indicates the exis- 
tance of a special condition in the data. A bit position in the 
" 0 " state carries no meaning other than that the data are'normal. 


The interpretation of the bit positions is as given below. 

Bit in H 1 ri Condition Implication 

1 Polarized power data questionable or 

no good 

2 Unpolarized power data questionable 

or no good 


3 

k 

5 

6 


Polarized noise level, n, (cf item 
37 ) changed this data record 

Polarized integration bounds B , B , 
(cf item 37 ) changed this data^rec 

Unpolarized noise level, n, (cf item 
39 ) changed this data record 

Unpolarized integration bounds, B , 
(cf item 39 ) changed this data 
record 


7 


System gain changed this data record 
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Explanation 
Bit Position "1” 

Polarized data may be flagged for any one of a number of reasons. 
The presence of interf erence, an error in setting the integration 
bounds, or a tape drive error in the data processing are examples 
of difficulties that would result in such a flag. In case of gross 
errors the flag represents an objectively known bad data point. 

In the case of a more subtle phenomenon suph as interference, the flag 
represents an experimenter* s subjective opinion. For the 116 cm data 
interference is the predominant cause of a data bad flags. Data 
users who wish to examine this question for themselves may do so by 
reprocessing the data from the JM Doptrack tapes. It is 
very strongly recommended that no flagged data be used without taking 
this precaution. 

Bit Position "2" 

Unpolarized data may be flagged for the same reasons as those 
given for the polarized data under bit position "1" above. However, 
because gross errors for polarized and unpolarized data may be inde- 
pendent, and because the unpolarized data possess a certain immunity 
to interference, which tends to be polarized, the flags in bit 
position "1" and " 2 " are not necessarily coincident. As before, 
in the case of interference, the investigator* s judgment is involved. 

Bit Position ,f 3 n 

This flag is set when the noise level, i.e., n in item 37, is 
changed during the data reduction process. This flag serves to 
alert the user that such a change has been made. Any discontinuity 
that occurs when this bit is set is likely to result from this cause. 
In the case of the polarized power such discontinuities are generally 
quite small, on the order of 1$. However, in the case of the 
normalized polarized power such discontinuities may be large, on the 
order of two. Obviously, no physical significance should be attached 
to such discontinuities. 
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Bit Position n 4 Tr 

Changes in the integration bounds are flagged for reasons simi- 
lar to those given under bit position " 3 "* The bounds, B^, B^, vary 
with the changing width and location of the polarized echo* Usually, 
such changes are very small and their effect is not noticable in the 
data* 

Bit Position rt 5" 

Changes in the unpolarized noise level, n, are made for the 
same reasons as those described under bit position ir 3 n * As before, 
changes in the unpolarized power, and especially the normalized 
unpolarized power, that occur with these changes in n are non-physi- 
cal * 

Bit Position TT 6 tT 

The comments that apply to the polarized integration bounds 
given under bit position ,r 4 M also apply here* 

Bit Position " 7 n 

System gain changes may occur during data reception in the 
receiving systems or during data playback in the record reproduc- 
tion system* Such changes are flagged since they will appear as a 
change in the polarized and unpolarized signal levels* Gain changes 
do not affect measures of echo bandwidths . 

48 * Antenna Gain is the gain of the spacecraft antenna in the O' * 8 

g 7 H'g 

direction determined from antenna patterns contained in NAA (1966b) 
(ll6 cm), MSC (1967) (Apollo ll+, 13 cm), or NAA (1969) (Apollos 15 and 
l6, 13 cm)* 

49 . Not Used 

5O* Sequence Number , where the first data record=l* 
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Cross-Reference Table for JM Doptrack/Integral Tape Trajectory Parameters 
The table below provides a cross-reference between the JM Doptrack 
trajectory parameter records and the Integral tape data records, 

JM Doptrack tapes are described in Appendix IV. That appendix gives the 
detailed format for those tapes, but does not define all of the trajec- 
tory parameters. Integral tapes are described in detail in Appendix VI. 
That appendix does include a detailed description of the various trajec- 
tory parameters. All of the trajectory parameters given on the JM Dop- 
track tapes are also found on the Integral tapes. This appendix provides 
a convenient means for determining the location of a given JM Doptrack 
tape trajectory parameter in the Integral tape format. For reference, 
the trajectory parameters are found in the sixth record of a data frame 
on the JM Doptrack tape (i.e., in record number N*6, N = 1, 2, 3, ..., 
where the first data record = l). 

JM Doptrack Integral Tape 

Word No. is found in Word No, 


1 

2 

3 

4 

5 

6 

7 

8 

9 


1 

2 
31 
30 

•28 

34 
9 

35 

36 


10 

n 

12 

13 


3 

4 

5 

13 


i4 


14 


15 

16 

17 

18 


15 

24 

25 
33 
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JM Doptrack Integral Tape 

Word No . is found in Word No . 

19 32 

20 26 


21 

22 

23 

2b 

25 

26 

27 


27 

29 

22 

23 

16 

17 

18 


28 

29 

30 

31 


6 

7 

8 

10 


32 


11 


33 


12 


Units and scale factors are identical on JM Doptrack and Integral 


tapes 
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Matrix Correction Factors 

The tables below give the matrix correction factors (c-matrix) that 
were applied to the ll 6 cm data in Subtask 5 (cf text). No corrections 
were applied to the 13 cm data* 

In all cases, the following values were used for C and C ^ 


C « 0.0 + j 0.0 

C = 1*0 + j 0.0 

22 


Only C„ „ and C „ were varied. Start times shown below refer to the 
11 21 

time on the first frame to which the particular c-matrix was applied; 
the same c-matrix was used on all succeeding frames with times less than 
the next entry in the table. 


Apollo lb 


Time UT2 (sec) 


start 



c n 

23838-35 

1.100 + j 

25155.63 

1.960 + j 

25166.28 

1.230 + j 

25168.94 

.1.100 + j 

25190.24 

1.200 + j 

25192.91 

1.960 + j 

25198.23 

1.200 + j 

25206.22 

1.100 + j 

26348.41 

0.929 + j 


C-matrix 


0.0 

0.0 

+ 

j 

0.0 

0.0 

0.050 


J 

0.200 

0.0 

0.070 

+ 

j 

0.150 

0.0 

0.0 

+ 

j 

0.0 

0.0 

0.0 

+ 

j 

0 

» 

0 

0.0 

0.050 

+ 

j 

0.200 

0.0 

0.0 

+ 

j 

0.070 

0.0 

0.0 

+ 

j 

0 

6 

0.0 

0.189 

+ 

j 

0.0 


Apollo 15 


all data 


C-matrix 


C 


11 


c 


21 


1.000 + j 0.0 


0.0 


+ j 


0.0 
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Apollo 16 


Time UT2 (sec) 




C-matrix 





start 

( 

-1 




r 





11 



21 


4592.35 

0.913 

+ 

j 

0.0 

-0.134 

- 

j 

0.045 

4780.77 

0.958 

+ 

J 

0.0 

-0.040 

- 

j 

0.010 

5432.35° 

0.944 

+ 

j 

0.0 

-0.055 

- 

j 

0.032 

5587.79 

1.000 

+ 

J 

0.0 

1 

0 

b 

0\ 

0 

+ 

j 

0.020 

5983.46 

0.900 

+ 

J 

0.0 

-0.100 

- 

j 

0.050 

6272.35° 

0.830 

+ 

J 

0.0 

-0.100 

- 

j 

o.o4o 

6484.32 

0.870 

+ 

J 

0.0 

-0.090 

+ 

j 

0.0 

7030.72 

0.860 

+ 

j 

0.0 

-0.060 

- 

j 

0 

lf\ 

0 

• 

0 

7H2.35 q ' 

0 

• 

CO 

f— 1 

+ 

3 

0.0 

-0.072 

- 

j 

0.490 

7333-74 

0.815 

+ 

3 

0.0 

-0.066 

+ 

j 

0.010 

7701.15 

0.778 

+ 

j 

0.0 

-0.060 

- 

j 

o.o4o 

7952.35° 

0.740 

+ 

j 

0.0 

-0.100 

- 

j 

o.o4o 


ct 

These changes are located in overlap areas described in Appendix IV (D3), 
The change takes effect on the frame with the time (word 2 , record 6; 
cf Appendix IV) closest to that shown here/ 
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Appendix IX 

Coordinate Transformations 

A* General Definitions 

r ss specular point position vector 

p = spacecraft position vector 

v = spacecraft velocity vector 

a „ = vector in local horizon system 

_AH 

a A , T = vector in vehicle reference frame 
_AV 

z= vector in primary or principal reference system 

0 ) 

Y / defined in C below (Euler angles between local horizon and 
0 / vehicle systems) 

look angles to specular point in vehicle system 
6 | orientation of vehicle along a y (3 direction (defined below) 

B. Data supplied by MSC gives spacecraft attitude in tT Local Horizontal” 
coordinate system. 

Definition of "Local Horizontal” system (see figure which follows) 



ZH 


YH 


XH 


~ P 
/P/ 

v x p 
/v X p/ 

x a ZH 
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ORBIT 


is directed toward the nadir, 

a YH iS directed a l° n g angular 
momentum vector, 

is in the plane of the orbit 
roughly along v. 


C. Specification of spacecraft attitude 

The Euler angles 0, Y ^ 0 specify the sequence of rotations required 
to transform Local Horizontal coordinates to Vehicle coordinates as 
follows : 



where the matrix 


ix | 

(axis ) J 


denotes a rotation of (ang) 
current (axis). Rotations are in the sequence (Y, Z, 
of rotation is clockwise when looking in the + (axis) 
instance^ consider: 


about the 

X) * The s ense 

direction. For 



6 7 



Appendix IX (cont.) 



= -X 2 sine* + X 3 cosa 


then the transformation is 



■0*i pi pi r* * 

.( X )J L(z)J L(y)J ” [ . 



cosY cos 0 sinY -cosY sin 0 

sin0 sine - cos0 sinY cos 0 cos0 cosY cos0 sinY sine + sin0 cos 0 

_sin0 sinY cose + cos0 sin 0 -sin0 cosY cos0 cose - sin0 sinY sine 


D. Conversion of Reference Coordinate to Local Horizontal System 

Quantities that are known in the principal reference system may be 
converted to the Local Horizon system through the directional cosines 
connecting the two. These are given by the vector components of the 


a AH, ex P ressed in the primary system. The a form the rows of the 
7 AH 
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transformation (rotation only) matrix. Thus 




a xn 1 


1 

on 


9 YH 2 

*yh 3 

>1 

a zn 2 

‘ ZH 3_ 



°YP 



Components in sub -H system are 
projections onto reference axes 
of that system. 
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E„ Conversion of Reference Coordinates for Vehicle Coordinates 

Conversion from principal coordinates to Vehicle coordinates may be 
made by two successive rotations, principal coordinates - — ►Local Horizon- 
tal coordinates Local Horizontal coordinates — ►Vehicle coordinates. 


or 


v] - M h - [ r vh] M [*p 


where R _ implies a rotation from the j to i system. Thus 


[v] = j*VH Vi. 


F. Computation of spacecraft to specular point look angles, a f p . 
1. The direction to the specular point from the spacecraft is 



In Vehicle coordinates 

'«] - [®vp] [v 

The quantity a is measured from the + X direction } 


a s = arccos(y xv ), 0 ^ Of ^ ir. 

The quantity is the azimuthal angle (minimum angle to plane) 

measured from the -Z^ axis, positive towards + Y^ r : 


P - signum 
s 

where signum (x)= X . 

I X I 

which follows : ' ' 


(Yyv) * arccos (“Yzv^ -Jt « P « JT, 

These quantities are depicted in the figure 
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P 



H . Sources 

Additional information regarding the definitions of Local Horizon 

and Vehicle coordinates may be found in MSC publications describing 

trajectory tapes ( MSC f 1970). Input for a , p f 6 computations are MSC 

s s 

supplied experimenter trajectory tapes. The principle reference frhme 
is selenographic . 
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Appendix X 

Relationship Between the Coherency Matrix and Other Specifications of 
Polarization 

The methods used to estimate the polarization state of Apollo 
bistatic-radar data are those described by Tyler and Simpson, ( 1970 ) . 
This appendix describes the relationship of the coherency matrix to 
the more common descriptors intensity, axial ratio and orientation of 
the polarization ellipse. 

The coherency matrix is from Born and Wolf , (1959). 

i x (t) = Re j E 1 e jtt)t J ^ 

e 2 (t) = Re j E 2 e JtUt | e £ 



J = 


<E^*> *> 

<E* E^> <E 2 E 2 > 


where are complex magnitudes associated 

with any pair of orthogonal linear polarizations. 


J 




J 


22 J 


J. 

ij 


<E.E* > 
1 3 
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= Tr ( j) 


~°11 

P 12 

- P 21 

P 22 - 


; Tr(j) - J 4- ^ no ~ received power 


Born and Wolf show that the percentage polarization y, i.e,, the fraction 
of Tr(j) that may be described by a deterministic polarization ellipse, 
is 

v = v'l-4 (p n 0 22 - P 12 o 21 )- 

Then any J may be written as 



unpolarized part polarized part 


where 

q . . = — (p . . - — ( 1-y )) ; q . . — 0 . . • 

11 y ' K 1 1 2 V T// i j y 1 J 

Ko (1962) gives the relationship between q and the Poincare sphere 
(Beckmann, I96S) as: (shown schematically in the figure which follows) 

q ^ = cos2o>, 0 £ ao ^ it 

j ?X 

q = sin^o cosco -) tan 2X 


1 m (q 12 ) 


_ * < 

* 2 


2 X < | 


Also note: 


tanto 



tanX > 0 


right elliptical polarization. 
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RCP 



7T/4 


For Apollo, a similar matrix where the sub— c designates decomposition 
of the incident wave into circular components, is used; 


<EE*> <EE*>" 

P A P. r 


J C11 J c 1 2 

where E and E 
P r 




denote left and 

<E E > <E E > 


J _ J 

“ 

l r r r J 


L c21 c 22 J 

right circular 


waves # 


Defining p and q 
-c -c 


in an analogous manner: 
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J -c = 


P cll P cl2 _ 1 
p c21 P c22 


(1-7) Tr(J c ) 


q q 

+ Y Tr(j ) c11 } where 

q c21 q c22- 


as before 


Y » - b (p ci 1 P 


c22 P cl2 P c21 


However, the interpretation of the q is modified in accordance with the 

c 

new definition: 


, v M ^ q cll + ^ q c22 
tanX - — = 

cn = | arg (q 1£ ), 

where 00 is referenced to instantaneous position at time t = 0, 

The invariance of Tr(*), Det ( • ) , and y under transformations to circular 
coordinates can be verified by direct computation. 
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XDS Sigma 5 Machine Images 

All data on the JM Doptrack and Integral tapes are in the form of 
XDS Sigma 5 machine images * One word is 32 bits, one byte, 8 bits. 

The data are either real (floating-point), integer, or alphanumeric. 

The machine images are as follows: 

A, Real 

A single precision real number consists of a sign bit (bit 0) , a 
biased^, base 16 exponent called a characteristic (bits l“7)j anc * 
a 24 bit fraction, A double precision number consists of a single 
precision number followed by an additional 32 bits of fractional 
significance (cf figure below). Unless otherwise noted, all numbers 
referred to in this report as "real" are single precision. 


+ 

Characteris- 

tic (c) 

/ 

1 

Fraction (F) 

— 1 

1 


0 

7 


f 

15 

23 

1 

31 

Extra Fraction Bits 
1 

1 

for Double Precision 

L 

1 

Only 

1 



32 63 


^The bias value of 64 jq is added to the exponent to make it possible 
to compare the absolute magnitude of two numbers without reference to 
a sign bit. , 
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A real number (n) is defined as follows: 

.064 

N = F Xl6 where 
F = 0 or 

2 ^|f|< 1 (single precision) or 

2 ^ <|f(< 1 (double precision) 
and 127 


A negative real number is the two 1 s complement of its positive 
representation, (Note that this differs from the IBM 360.) 


B . Integer 

Integers are written in full-word, two 1 s complement representation. 


C, Alphanumeric 

Alphanumeric data, strings of characters, are represented in Extended 
Binary-Coded-Decimal Interchange Code (EBCDIC), Each character 
occupies one byte (8 bits), so each word of alphanumeric data con- 
tains four characters. 
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NOTES 


Block Diagram I: Block Diagram of 

Stanford Apollo Bistatic-Radar Data Processing and Reduction 


(1) a) Receiving system described in Appendix I * 

b) Critical filter response given in Appendix III. 

c) RCP, LCP, coherency maintained through entire system. 

d) Antenna pointing maintained toward center of moon using 
Stanford Research Institute lunar ephemeris. 

e) Absolute phase in LCP, RCP channels not controlled, relative 
phase maintained. 

f) Spectral purity of downlink signals measured for Apollo 
equipment type ( MSC , 1970 ). 

(2) a) Deep space network station under NASA control, used standar- 

dized procedures except as noted, 

b) Apollo operational system used for 64 m dish pointing, closed 
loop receiver acquisition. Station configuration for bistatic- 
radar experiment described elsewhere ( DSN , 1970, 1971). 

c) Absolute phase in RCP, LCP not controlled, relative phase 
maintained, 

( 3 ) a) Stanford Signal Conditioning Unit provided critical control 

over system frequency response. Filter characteristics 
given in Appendix III. 

b) FR1400 tape recorders A, B used to provide continuous data 
across tape changes. 


(4) a) H.P. 3955 B frequency response adjusted for a maximally flat 
frequency response on playback using calibration tapes from 
data source machine (H.P. 3955 B for ll6 cm data; FR1400 A, B 
for 13 cm data) . 

b) Coherent sampling maintained. Reference signal from tape 
used to synchronize samples. 


c) Actual time recovered by use of time code translator clock 
output to establish start times. Time from start main- 
tained by counting sampling pulses. 
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(4) d) 


Actual sampling rates: (cont.) 


Ratio of Effective 
Sampling Rate to 


Flight 

Wavelength 

Sampling Rate 

Actual Sampling Rate 

Apollo 14 

13 cm 

10.75 

kHz 

4 

Apollo 14 

ll6 cm 

10.00 

kHz 

1 

Apollo 15 

13 cm 

10.75 

kHz 

4 

Apollo 15 

ll6 cm 

10.00 

kHz 

1 

Apollo 16 

13 cm 

21.5 

kHz 

2 

Apollo l6 

ll6 cm 

10.00 

kHz 

1 


( 5 ) a) Fourier coefficients computed from successive groups of weighted 
data: data group length either 1024 or 2048 samples, weighting 

function is sin^ (jtt/T)^, where t is time and T is duration of 
sample group length . 

b) Data analysis lengths 

Data Source Sample Length Analysis bandwidth (Hz) 


Apollo 14 S-band 

1024 

42.0 

Apollo l4 VHF 

1024 

9.8 

Apollo 15 S-band 

1024 

42.0 

Apollo 15 VHF 

2048 

4.9 

Apollo l6 S-band 

1024 

42.0 

Apollo l6 VHF 

2048 

4.9 


(6) Polarimeter computed elements of signal covariance matrix 
(see Subtask 5) . 

(7) a) Normalizing data obtained from data runs employing noise input. 

Purpose was to compensate for variations in receiver passband, 

b) Normalizing tapes contain smoothed power spectra from noise 

source data. 

(8) a) C^matrix correction was polarization coordinate transformation 

to correct for errors in receiver antenna system. Applied to 
ll6 cm data only. 

b) Criterion for choosing c-matrix was to minimize polarized part 
of system noise. 
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( 9 ) a) Merge combined data tapes with geometrical parameters obtained 
from MSC (Best-Estimate-Trajectory) Post-Plight tapes. 

b) MSC data are interpolated to center of data averaging window. 

c) JM Doptrack tapes constitute primary data source for analysis. 

(10) Carrier suppression deleted direct signal from the data by 
means of an empirical algorithm (see Subtask 7 and Appendix v). 

(11) a) Tapes edited for changes indicated in operational logs and 

interference. 

b) Handscaled data added to tape from card source. 

c) Edited integral tapes constitute reduced data records. 
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